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I.  DESCRIPTION  OF  PROBLEM  AND  RESEARCH  OBJECTIVES 

The  ability  of  atomic  gases,  or  polar  molecular  gases  oriented  by 
electric  fields,  to  double  and  triple  the  frequency  of  laser  radiation  is 
determined  by  the  hyperpolarizability  tensors  B  and  Y  of  the  gaseous  mole¬ 
cules.^  Thus,  these  tensors  are  becoming  increasingly  important  in  a  wide 

12-41 

variety  of  real  and  potential  laser  applications.  '  Such  applications  in¬ 
clude  the  possibility  of  producing  high  frequency,  high  power  laser  radiation, 

(2-5) 

utilizing  various  up-conversion  techniques. 

For  example,  it  has  been  proposed  to  use  selenium  vapor  to  fre- 

>  (4) 

quency  triple  CO2  laser  radiation  to  enable  atmospheric  transmission'  and 

it  has  been  shown  that  up-conversion  of  infra-red  laser  radiation  utilizing 

(2) 

two-photon  pumped  alkali-metal  vapor  is  possible.  '  Since  the  crucial  piece 
of  information  in  each  case  is  the  hyperpolarizability  of  the  vapor,  the  prac¬ 
tical  development  of  such  techniques  is  jeopardized  by  a  lack  of  knowledge  of 
the  hyperpolarizabilities  of  atoms  and  molecules.  Experimental  determination 
of  these  quantities  are  difficult,  and  the  few  theoretical  calculations  which 
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had  been  made  prior  to  the  current  project,  were  not  able  to  even  provide 

order  of  magnitude  estimates  of  these  properties. 

Hyperpolarizabilities  may  be  deduced  experimentally  by  direct 

observation  of  harmonic  generation  in  gases^’®’^  and  from  measurements  of 
(8  9 1 

the  Kerr  effect.'  ’  '  The  experiments  are  complicated  and  the  range  of  un¬ 
certainty  is  sometimes  large.  Since  hyperpolarizability  is  an  atomic  or 
molecular  property,  it  may  also  be  predicted  from  quantum  mechanical  considera¬ 
tions.  Thus,  it  is  quite  important  to  have  reliable  theoretical  predictions  of 
hyperpolarizabilities  to  complement  the  experimental  effort,  and  to  provide 
results,  particularly  for  hyperpolarizability  as  a  function  of  frequency,  for 


experimentally  inaccessible  problems.  _ 

*^he  objectives  of  this  effort  fer-AF&Sf/ are  threefold: 

I ^  Develop  state-of-the-art  quantum  mechanical  methods  to  predict, 

from  first  principles,  hyperpolarizabilities  of  atoms  and 


moleculesj 

Ilj  Determine  why  previous  theoretical  calculations  have  not  been 
able  to  even  provide  order  of  magnitude  estimates  of  experi¬ 
mental  hyperpolarizabilities^ 

^III^  Apply  the  quantum  mechanical  techniques  developed  in  this  work 
to  predict  hyperpolarizabilities  and  non-linear  susceptibilities 


for  a  number  of  atoms  and  molecules  of  interest. 


This  effort  will  provide  complementary  information  to  on-going  experimental 
efforts  directed  toward  obtaining  hyperpolarizabilities. /L 

To  improve  the  theoretical  predictions,  four  predominant  areas  have 
been  identified.  These  are:  (1)  the  basis  set  problem;  (2)  the  effect  of 
correlation;  (3)  the  sensitivity  of  the  hyperpolarizabilities  to  geometric 
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changes  due  to  vibration;  and  (4)  the  frequency  dependence.  In  this  study  we 
have  mostly  concentrated  on  investigating  (1),  (2),  and  (3)  with  the  question 
of  frequency  dependence  to  follow  once  the  initial  three  questions  have  been 
resolved.  Once  an  understanding  of  the  elements  needed  to  make  accurate 
calculations  is  achieved,  then  a  series  of  large-scale  computational  studies 
of  several  atoms  and  molecules  of  special  interest  could  be  undertaken. 

In  our  initial  study,  we  have  developed  computer  programs  needed  to 
obtain  correlated  hyperpolarizabilities.  ^  We  have  used  these  programs  to 

(10-13' 

provide  the  first  correlated  calculations  of  hyperpolarizabilities  for  molecules, 

we  have  made  the  first  study  of  the  dependence  of  the  hyperpolarizability  on 
(12) 

molecular  geometry;  and  we  have  determined  the  type  of  basis  sets  that  are 

needed  to  provide  reliable  results.  ^ 

Armed  with  these  components,  our  initial  calculations  for  the 

hyperpolarizabilities  for  molecules^10-13^  provide  significantly  better  agreement 

with  the  observed  experimental  value^^than  any  other  theoretical  work  has  been 

able  to  achieve.  It  has  been  clearly  demonstrated  that  correlation  effects  are 

of  enormous  importance  in  ab  initio  theoretical  predictions  of  hyperpolarizabilities. 

For  example,  even  in  the  comparatively  simple  case  of  NH^  correlation  changes 

8  by  a  factor  of  four,  and  doubles  the  sum  8=(e  +8  +6  ).  In  more 

£££  lll  AAt  yyy 

polar  molecules  such  as  BeO,  the  correlation  effect  can  be  even  more  dramatic. 

Equally  large  effects  are  found  in  v  .  Also,  the  extreme  geometry  dependence 
of  j3  and  ^  has  been  elucidated  for  the  first  time.  These  results,  plus  a  number 
of  additional  accomplishments  in  the  program,  have  set  the  stage  for  a  rigorous, 
predictive  theoretical  approach  for  hyperpolarizabilities,  that  should  serve 
to  assist  in  the  understanding,  and  in  the  experimental  interpretations  of 
nonlinear  optical  processes  in  gases. 
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II.  HYPERPOLARIZABILITIES  AND  NONLINEAR  OPTICS 

The  treatment  of  optical  effects  in  a  molecule  is  based  upon 

the  dipole  moment  induced  in  the  medium  when  perturbed  by  monochromatic 
(151 

radiation.  The  appropriate  harmonic  perturbation  operator  is 


H^(r,t)  =  -  ej  r^E^cos  ct  ,  (1) 

=  xh(1)(r)(eiwt+e_1'ut)  ,  (2) 

subject  to  the  definitions 

X  =  |E(w)!/2  (3) 

(D  N 

h*  (r)  =  -e  l  r}  a  .  (4) 

1-1  “ 


In  Equations  (1)  and  (4)  N  is  the  number  of  electrons,  E/^the  field 
strength  vector  for  frequency  u,  r^  the  position  vector  for  each  electron  i, 
and,  a,  a  unit  vector  in  the  polarization  direction. 

The  induced  dipole  moment  is  defined  as 

jj.(t)  =  e<f|r|’?>  (5) 

where  f(r,t)  are  the  eigenfunctions  of  the  time-dependent  Schrodinger  equation. 
Normal  considerations  of  time-dependent  perturbation  theory^  provide  the 
expansion 

f(r,t)  ■  f^(r,t)  +  f^(r,t)  +  ’?^(r,t)  +  ... 


(6) 
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By  separating  the  spatial  and  time  dependent  parts  of  the  perturbed 
wavefunctions,  we  may  obtain  from  Equation  (5), (15) 


<w|rj4->=P(0)+P(1)EU)cos  „t+P^2)(E(w))2"P(2)(EU))2  cos  2  ut 

1  — 1  —0  — fjj  — u  — 

+  p(3)(r(“))3COs  ut  +  P(3)(E(w))3  COS  3  ut 


(7) 


The  coefficients  in  Eqn.  (7)  are  electric  susceptibilities.  The 
terms  that  are  quadratic  and  higher  in  the  field  strength,  are  responsible 
for  the  non-linear  optical  effects  of  a  medium,  as  indicated  by  the  frequency 
doubling  and  tripling  terms,  cos  2  ut  and  cos  3  ut.  The  main  objective  of 
the  present  study  is  to  investigate  the  capability  of  first-principle  quantum 
mechanical  predictions  of  these  susceptibilities  for  molecular  gases. 

The  definitions  of  the  susceptibilities,  P^n\  jn  terms  of  the  SDatial 

mu 

functions  are 


(8) 


Pf,  ^  +«/^r|<|> 


(°). 


l/4[<^°)!r|^2)>  +  <^2)|r|^°)>  +  <4i1)|r|4J)> 


p|2)=  l/2[«>o(0) \r\t[ZJ>  +  <^|r|*J0)>  +  Ini ♦i^> 


(9) 

(10) 


(ID 
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-<*l0)ir,*«oui;>i4;>>3  . 

P«3).  l/4C<*i°)!rUii3>,^^r|^  *<♦«;>  M<}> 

+42 }  iii  ♦!  | }> + <^2 }  iii  + <*! j  *  iri  42  }> 

-  <4° }  i  m  4° }  >  <  <4; 5 1 4: } }  I  ^ 5 1 4; }  > 
^>j2,>).<>irurw  )in/1i)>K<^;)i4;)= 

+^,i<w;>u!;)>)]  . 


^,|r|*i!)»-«*‘,,M*i,)*(<^*^> 

^;)i<)»-</;)i<)>«*'o)iti*<;,>^;>iii*io)>) 


(15) 

These  terms  have  convenient  physical  interpretations.  Equation  (3) 
is  the  static  dipole  moment  while  Equation  (9)  is  the  frequency  dependent 
dipole  polarizability,  which,  of  course,  is  also  related  to  the  index  of  re¬ 
fraction,  The  remaining  terms  are  qenerally  referred  to  as  hyperpolarizabilities, 

(2) 

or  nonlinear  susceptibilities.  P^  ,  the  first  nonlinear  term  in  the  expansion 

12) 

of  the  induced  dipole  moment,  determines  optical  rectification.  P'  defined 

— 2<a) 

in  Equation  (11),  leads  to  second-harmonic  generation.  The  two  third-order 
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terms  are  respectively  responsible  for  the  intensity  dependent  refractive 
index  and  third-harmonic  generation.  By  generalizing  the  perturbation  to 
include  a  combination  of  ac  and  dc  electric  fields,  additional  phenomena 
like  the  linear  electro-optical  effect,  dc-induced  second-harmonic  generation, 
and  the  Kerr  effect  may  also  be  described. 

Equations  (8-1 3)  are  the  most  appropriate  expressions  to  employ 
in  actual  calculations  of  the  coefficients  in  Equation  (7),  since  the 
wavefunctions  are  defined  by  variational  equations  obtained  from  time- 
dependent  perturbation  theory.  However,  for  additional  insight,  particularly 
with  regard  to  static  (ui=o)  hyperpolari2abilities,  it  is  often  easier  to 
think  in  terms  of  sum  over  states  perturbation  formulae P^'^Hiith  tfie  definition 
of  h^(r)*h(r)  in  Equation  (4),  and  ukQ  being  the  excitation  energy  to  state  k. 
Equations  (9-11)  become. 


■  *J0  !WVo/t(“kot"W“ko-“n  ' 


(14) 


sl2)  -  -1/4  I(2rok\0h00/“ko['“kot“,:l  +(“ko-“)>  Ioohokhko/t(“ko'f“)2*(,-ko"“)2:,! 


-0 


k^o 


+1/4k,i^0  {2rokh^o/ukoC(^o+“)+(^o'U,):! 


(15) 


+hok^£heo/C(“ko+u)(^o+^+(“ko-“)(^o-w)] }  ’ 


=  -1/2.1  <^if^0^n/C(“l(o+2u^UkO+")  +  (l0ko-2u)Uko'u);l 


— 


Wo 


06) 


*-Coohokhko/[(“ko+“)  (“ko-“)]) 
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oo 


+V2k  ^0{^kh^VCUko+2“)U£o+''°)+Uko'2u))(u£c^' 


+  hokrk£h£o/^u3ko+u)^a,£o',J)^ 1 

The  expressions  for  the  third-order  coefficients  are  very  long  and  will  not  be 
reproduced  here  [see  Reference  (15)  for  these]. 

Another  notation  is  frequently  used  for  the  susceptibilities  in 
Eqn.  (7),  particularly  for  the  static  U=o)  case.^  The  defining  expression 
is  a  Taylor's  series  expansion  of  the  energy  in  the  presence  of  an  external 
field, 

u(|)  -  M(o)  -  UjE,  -  (2! )*'  (17) 

-  <3!>_1  Mjk  E(  Ej  Ek 

-  I41*"'  ’SjW  E,  Ej  Et  Ef 


The  summation  over  repeated  indices  is  assumed.  The  quantitives  u-j ,  a^, 
B-jjk,  and  Yijk£  are  respectively  components  of  the  permanent  dioole  moment, 
polarizability,  and  the  hyperpolarizabilities.  A  derivative  of  Eqn.  (17) 
with  respect  to  a  component  of  E.  defines  the  expansion  of  the  induced  dipole 
moment  which  may  be  compared  to  Eqn.  (7).  Some  consideration  of  Eqn.  (8) 
and  Eqns.  (14-16)  shows  that  in  terms  of  sum-over-state  expressions, 


(20) 


Eqn.  (20)  and  similar  expressions  for  are  used  to  make  some  comparisons 

with  the  frequency  dependent  results  for  P_^  and  £.3^ •  Since  tne  excitinc 
frequency  u  is  expected  to  be  comparatively  low,  often  in  the  infra-red,  or  visible 
static  predictions  for  6  and  tr  are  not  expected  to  differ  substantially  from 
the  frequency  dependent  experimental  results.  However,  by  considering  an 
expansion  about  u  in  Eqn.  (16)  a  more  formal  expression  for  the  frequency 
dependence  can  be  obtained. 

Using  the  relationship  that 

(a+b)"^  =  a"^  t  a~^ba"^  +  a'^ba^ba”1  t  .  .  .  (21) 

to  terms  quadratic  in 


£<1!  =  1/23  +  (1/2  l 


k,f. 


ok-Mc?  to  ko  to' 


ko 


'to' 


(22) 


+  7  M  rokhkth£o/a,kow;,o(“/u,k0)(w/^o)} 


Eqn.  (22)  demonstrates  that  all  terms  linear  in  (w/ukQ)  vanish.  For  a 
frequency  to,  chosen  in  the  visible,  the  ratio  of  |<*>/wk  |»  when  u|<o  in  an 
optical  transition,  is  usually  small.  Hence,  the  quadratic  dependence  also 
attests  to  the  comparatively  small  difference  that  would  be  expected  between 


>(2)  a  nr!  1/9(3 
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Defining  ^  as  the  lowest  optical  excitation  energy  in  the  molecule, 
then  |u)/u)|2|a)/iJko|  for  any  k,  and  |u/Q2  can  be  factored  from  Eqn.  (22).  If 
the  two  contributions  to  8.  have  the  same  or  different  signs,  then  it  also 
follows  that  [ 7 s ( ui/aT)  j  2  l  the  second  term  in  Eqn.  (22).  Hence,  some  estimate 
of  the  frequency  dependence  in  6_  can  be  made  using  this  expression.  This  gives 
a  percent  error  in  the  static  s_  values  to  be  ^  7(w/uT)2  x  100. 

The  main  thrust  of  this  research  program  is  directed  at  ab  initio 
quantum  mechanical  predictions  of  hyperpolari ^abilities.  These  quantities 
are  difficult  to  obtain  experimentally,  and  equally  difficult  to  compute  from 
first  principle  quantum  mechanics.  The  next  section  discusses  the  quantum 
mechanical  approach. 
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III. 


QUANTUM  MECHANICAL  PREDICTIONS  OF  HYPERPOLARIZABILITIES 


A.  Computational  Approach 


The  quantum  mechanical  approach  to  calculating  static  hyperpolariz¬ 
abilities  is  built  upon  Eq.  (17).  Either  of  two  routes  to  the  evaluation  of 
the  coefficients,  ory.^  could  be  taken.  One  approach  is  to 

evaluate  W(E_)  at  a  series  of  small  values  of  the  field  strength  (i.e.,  the 
finite-field  method),  from  which  a  series  of  finite  difference  equations  can 
be  used  to  obtain  the  dipole  moment  and  polarizabilities.  The  other  approach 
would  directly  calculate  the  coefficients  by  considering  expressions  obtained  from 
the  appropriate  order  of  perturbation  theory.  Althouqh  the  second  route  does  not 
require  taking  small  differences  among  large  numbers,  it  is  more  difficult  to 
implement  at  a  high  level  of  sophistication  since  separate  sets  of  involved 
computer  programs  must  be  developed  for  each  order  of  perturbation  theory. 

The  finite-field  method  instead  supplies  all  the  coefficients  from  a  single 
computer  program  by  simply  evaluating  W(£)  at  enough  field  strengths.  This 
also  offers  the  bonus  that  a  much  higher  degree  of  electron  correlation  can 
be  incorporated  into  the  computer  programs  that  evaluate  W(E)  than  would  be 
convenient  to  develop  if  each  quantity,  u ,  £,  and  ^  were  to  be  evaluated 

separately. 

In  the  finite-field  approach,  there  is  also  a  possibility  that  a 
lower  order  coefficient  can  be  contaminated  to  some  degree  by  the  higher 
coefficients.  To  reduce  this  possibility,  one  can  choose  the  finite-field 
strengths  symmetrically,  such  as  +0.01  so  that  odd  and  even  combinations  of 
W(+0.01)  can  be  made  to  eliminate  the  next  term  in  the  expansion  from  the 
finite  difference  equations.  This  eliminates  most  of  the  possible  contamination. 
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The  finite-difference  equations  derived  for  u,a,£,  and  y  are  reported  in 
Reference  (10),  included  as  Appendix  A  to  this  report. 

The  energies  W(EJ  can  be  obtained  from  the  linked  diagram  theorem 
of  MBPT. (18"22>  This  is 


W(E)  =  Wchf(E)+£  <$0(I)|V(E)[(Wo(E)-F(E))'1V(E)]K|$0(ij>  (23) 

k=l  1 


when  the  subscript,  L,  indicates  exclusion  to  1 inked-diagrams.  The  terms 
in  this  expression  are  WCHp(Ej,  which  is  the  coupled  Hartree-Fock  finite-field 
value,  $Q(j[),  the  CHF  wavefunction,  and  F(E)  the  sum  of  one-electron  Fock 
hamiltonians  including  the  external  one-electron  perturbation  of  Eq.  (1). 

Wq(£)  is  the  sum  of  the  field  dependent  one-electron  energies  associated  with 
the  orbitals  composing  $0(E).  It  is  apparent  from  Eq.  (23)  that  the  zeroth- 
order  energy  in  this  expression  is  the  CHF  result.  In  previous  work  on 

(23  24) 

hyperpolarizabilities  of  molecules,  seldom  has  even  the  CHF  model  been  used;  ’  ' 

instead  various  uncoupled  approaches  (UHF)  have  usually  been  applied^25)  due 
to  the  computational  difficulties  encountered.  The  present  predictions  start 
with  CHF  and  include  the  predominant  correlation  corrections  arising  from  the 
linked-diagram  expansion  in  Eq.  (23).  In  this  respect,  this  work  is  unique 
in  studies  of  molecular  hyperpolarizabilities. 

In  addition  to  ignoring  electron  correlation  effects  in  previous 


predictions  of  hyperpolarizabilities,  the  basis  set  problem  has  also  received 
inadequate  attention.  In  a  molecular  calculation,  the  basis  set  is  usually 
chosen  to  be  composed  of  functions  (called  Slater  orbitals)  that  are  similar 
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to  atomic  orbitals  of  the  general  form  rn’^e*arY  Je,  $),  where  r,  0,  and  are 
the  usual  spherical  coordinates,  and  n,  1,  and  m  are  the  normal  quantum  numbers. 

The  exponential  scale  factor,  a,  referred  to  as  the  orbital  exponent,  reflects 
the  screening  of  an  electron  in  a  given  orbital  by  the  interior  electrons  in 
the  atom  or  molecule.  In  practice,  when  studying  polyatomic  molecules,  an 

_  gp2 

exponential  dependence  of  the  form  e  is  usually  preferred  because  of 

dramatic  simplifications  in  the  computation  of  molecular  integrals.  When  such 

functions  are  used  to  represent  a  Slater  type  function,  they  are  referred  to  as 

a  contracted  Gaussian  type  orbital  (CGTO).  In  this  work,  the  latter  functions 

are  used,  with  a  designation  for  NH^  perhaps  of  (5s4p3d|4s2p) meaning  5s,  4p,  and 

3d  type  orbitals  are  used  on  the  heavy,  nitrogen,  center,  and  4s  and  2p  functions 

on  each  hydrogen  atom.  All  components  of  a  p-function  (i.e.,  p  ,  p  ,  p  )  or 

x  y  z 

a  d-function  (d  ,  d  ,  d  ,  d  ,  d  ,  d  }  are  included  when  a  single  p  or  d 

jrjr  Ll.  Ajr  j L.  Aj 

is  listed  in  the  designation.  Hence,  {5s4p3d|4s2p) corresponds  to  65  CGTO 
functions  for  NH3. 

Basis  sets  have  evolved  over  the  years  in  molecular  quantum  mechanics 
until  sets  of  orbital  exponents  for  the  different  atoms  composing  the  molecule 
have  become  established  '  ;  for  ordinary  energy  related  properties.  That  is, 

the  principal  term  in  a  field-free  hamiltonian  is  proportional  to  1/r,  so  these 
basis  sets  reflect  that  dependence.  When  dipole  moments,  polarizabilities  and 
hyperpolarizabilities  are  wanted,  the  operator  dependence  is  instead  r11,  where 
n  *  1,2,  or  3.  This  places  a  stringent  requirement  on  the  basis  sets  to  be  used. 
Since  the  molecule  itself  must  be  adequately  described,  a  normal  (1/r)  type 
energy  basis  is  needed,  while,  in  addition,  a  description  of  the  long-range 
dependence  of  the  charge  density  is  crucial  to  accurately  represent  an  rn  operator. 
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Since  quantum  mechanical  applications  have,  until  recently,  been  rather 
limited  in  the  number  of  basis  functions  that  could  be  used  (typically 
<_50  CGTO's  at  the  CHF  level,  and  less  if  correlation  is  needed),  the  size 
of  basis  which  would  have  the  combination  of  properties  required  to  adequately 
describe  hyperpolarizabil ities  has  been  difficult  to  achieve. 

In  order  to  have  reasonable  confidence  in  the  capability  of  a 
basis  set  to  simultaneously  describe  the  electron  density  in  the  energy 
rich  (i.e.,  1/r)  region  of  the  molecule,  and  in  the  tail  (i.e.,  rn)  region, 
our  initial  molecular  calculations  studied  the  prediction  of  hyperpolari¬ 
zabilities  with  a  series  of  basis  sets  for  HF  by  comparing  with  numerical 

(27) 

Hartree-Fock  calculations.  '  The  latter  calculations,  which  are  only 
possible  for  the  u2,  azz>  6ZZZ>  and  Yzzzz  components,  avoids  any  basis  set 
error  since  the  differential  equations  of  Hartree-Fock  theory  are  solved 
numerically.  Besides  being  limited  to  the  z-components  of  the  various 
quantities,  this  technique  is  applicable  to  atoms  or  diatomic  molecules,  so 
its  value  lies  in  selected  comparisons  rather  than  general  use.  In  this 
manner,  we  are  able  to  show  that  a  type  of  selection  scheme  for  CGTO  basis 
sets^’^  appears  to  be  capable  of  reducing  the  basis  set  error  for  CHF 
hyperpolarizability  predictions  to  ^10  percent.  Basis  sets  so  chosen,  however, 
require  typically  65-80  functions  even  for  the  small  molecule  isoelectronic 
sequence  HF,  HgO,  and  NHj.  Such  large  numbers  of  basis  functions,  still 
sorely  test  quantum  mechanical  methods  and  the  efficiency  of  computer  programs, 
so  it  is  absolutely  necessary  that  the  more  sophisticated  and  computer 
intensive  correlated  calculations  benefit  from  the  most  advanced  methods  possible. 


The  development  of  MBPT^‘22)  under  AFCSR  sponsorship, HO, 28-30)  has  mac)e 
such  correlated  calculations  possible.  (The  reader  is  referred  to  Appendices  A 
and  B,  for  more  detailed  information  about  basis  set  selection  and  MBPT 
methods  for  evaluating  correlation  corrections  for  hyperpolarizabilities.) 

Two  other  theoretical  developments  in  this  research  should  also  be 
mentioned.  These  developments  pertain  to  infinite-order  correlated  calcula¬ 
tions  and  the  frequency-dependent  problem. 

Most  of  our  computations  in  the  two  years  of  this  program  were 

limited  to  SDQ-MBPT(4).  The  acronym  means  the  inclusion  of  all  correlation 

effects  due  to  single-,  double-,  and  quadruple-excitation  diagrams  that  arise 

through  fourth-order  in  the  linked-diagram  expansion.  This  model  is  described 

in  detail  in  Appendix  A,  and  we  have  established  this  model  to  be  a  highly 

(it) 

accurate  approach  for  ab  initio  correlated  calculations.  The  sensitivity 
of  hyperpolarizabilities  to  correlation  effects,  however,  suggested  that 
terms  higher  than  fourth-order  might  be  important  for  such  sensitive  properties. 

To  address  this  question,  a  substantial  development  in  the  past 
year  has  been  to  generalize  the  theory  and  computer  programs  to  carry  out 
CCSD  calculations  (i.e.,  coupled-cluster  singles  and  doubles).  This  model 
is  an  infinite  order  generalization  of  SDQ-MBPT(4),  to  which  CCSD  reduces 
in  fourth-order,  permits  a  consideration  of  possible  contributions  of  higher 
order  terms.  The  computer  codes  are  functioning  and  we  have  been  making 
applications  at  this  level  for  the  HgO  molecule,  but  the  excessive  computer 
time  requirements  have  prohibited  us  from  obtaining  a  complete  picture  at 
this  time.  Our  current  indication  is  that  little  change  will  occur  due  to 
these  higher-order  corrections,  but  we  will  reserve  judgement  until  our 
study  is  completed. 
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The  other  theoretical  development  we  have  been  pursuing  pertains 
to  determining  hyperpolarizabilities  as  a  function  of  frequency.  This  is 
being  accomplished  by  employing  time-dependent  Hartree-Fock  theory  (also 
known  as  the  Random  Phase  Approximation)  and  calculating  the  8  hyperpolari¬ 
zabilities  at  a  range  of  frequencies.  This  approach  has  the  bonus  that  a 
great  deal  of  information  about  the  excited  state  of  atoms  and  molecules  is 
also  obtained  from  inspecting  the  resonant  frequencies.  So  far,  we  have 
assumed  that  the  correlation  corrections  to  the  frequency  dependent  polari¬ 
zabilities  will  be  essentially  the  same  as  in  the  static  case.  Otherwise, 
the  computations  would  become  intractable.  We  have  made  some  preliminary 
applications  at  this  level  which  so  far  support  the  assertion  that  the 
frequency  dependent  hyperpolarizability  is  not  too  different  from  the  static 
value  unless  one  is  near  a  reconance  when  more  accurate  computational  approaches 
would  be  required.  Like  the  infinite-order  CCSD  approach  discussed  above, 
this  project  has  necessarily  been  aborted  by  AFOSR's  untimely  failure  to 
continue  this  research  project. 
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B.  Applications 

The  molecules  studied  in  the  two  years  of  this  project  include  Be, 

HF,  H2O,  NHg,  BeO,  LiF,  and  Ng.  In  Appendices  A  and  B,  HF  and  are  considered 
in  depth.  and  Be  are  discussed  in  previous  reports,  so  the  present  discus¬ 

sion  will  primarily  focus  on  NH^  with  a  few  additional  comments  about  comparisons 
among  different  molecules. 

From  the  viewpoint  of  a  comparison  between  experiment  and  theory, 

NH3  is  a  very  interesting  system.  From  dc-induced  second-harmonic  generation 
experiments.  Miller  and  Ward  find  the  electric  susceptibility  x/-|2^=  10/3(3ZZZ 
+  3  +6  )  to  be  -209+5  xio"33  esu/molecule.  The  previous  theoretical 

4aa 

calculations  at  the  CHF  level,  predict  -44.4 ^-65.1 ^-19.0^4^-15.6^23^  and 
(23) 

-40.8)  all  is  very  poor  agreement  with  the  experimental  result.  Besides 
the  poor  agreement  with  experiment,  as  a  basis  set  improved,  the  theoretical 
calculations  show  a  tendency  to  be  in  progressively  worse  agreement  with 
experiment.  This  is  indicated  by  the  fact  that  the  CHF  calculations  of 
Lazzaretti  and  Zanasi  obtained  -65.1  with  a  simple  double-zeta  basis  of  just 
16  CGTO's,  while  getting  -19.0  when  they  include  polarization  functions  on  all 
atoms  which  provides  a  30  CGTO  basis  set.^^  This  behavior  is  a  common  trait 
among  hyperpolarizability  calculations,  reflecting  the  fact  that  a  small  basis 
set  has  a  sufficiently  poor  description  of  the  electron  density,  that  gross 
error  cancellation  can  occur,  with  the  errors  in  CHF  theory  partially  opposing 
the  errors  in  the  basis  set.  Once  the  basis  set  is  somewhat  improved,  so  that 
a  better  discretion  of  the  density  is  possible,  the  deficiencies  in  the  CHF 
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model  become  more  apparent,  resulting  in  even  poorer  agreement  with  experiment. 
However,  a  30  CGTO  basis  is  still  in  no  sense  adequate  for  the  calculations 
of  properties  as  sensitive  as  hyperpolarizabilities.  In  our  work,  much  better 
basis  sets  are  used  while  electron  correlation  is  also  properly  included. 

Using  the  procedure  described  elsewhere^)  which  enabled  us  to 
define  basis  sets  that  provide  good  agreement  with  numerical  Hartree-Fock 
calculations,  we  have  generated  a  (5s4pld|4s2p)  basis  set  for  NH^  (basis  A; 

53  CGTO)  and  a  (6s4p3d |4s3p)  basis  set  75  CGTO.  The  results  obtained  with 
these  basis  sets  are  shown  in  Table  1.  A  susceptibility  of  "209  XlO  33 
esu/molecule  corresponds  to  a  8  =  -72  in  atomic  units. 

It  is  apparent  from  Table  1  and  II  that  the  much  superior  basis  sets  of 
the  present  work  are  capable  of  quite  accurate  predictions  of  the  dipole 
moment  and  polarizability.  However,  even  for  these  properties,  which  are  not 
as  sensitive  to  correlation  as  are  the  hyperpolarizabilities,  there  is  a 
significant  change  between  the  noncorrelated  (CHF)  results  and  the  correlated 
SDQ-MBPT(4)  predictions.  The  dipole  moment  changes  about  6  percent,  just  as 
in  HF  and  H2O,  while  the  polarizability  changes  by  ^9  percent.  For  these 
two  properties,  there  is  far  less  difference  between  the  two  basis  sets  than 
there  is  due  to  correlation  effects.  Notice  in  particular  in  Table  II  the  dramati 
change  in  the  anisotropy,  due  to  correlation. 

The  most  dramatic  effect  of  basis  set  and  correlation  is  reserved 
for  the  more  sensitive  hyperpolarizabilities.  For  example,  the  change  to  a 
more  complete,  diffuse  basis  set,  causes  a  factor  of  two  reduction  in  ezzz 
which  is  responsible  for  basis  3  predicting  a  |s|  that  is  smaller  than  in  the 
somewhat  poorer  basis  A,  but  still  the  agreement  between  the  CHF  results  and 

,'i 
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TABLE  I.  COMPARISON  OF  CHF  AND  SDQ-MBPT(4)  DIPOLE  MOMENTS 
AND  POLARIZABILITIES  FOR  MOLECULES  (ATOMIC  UNITS) 


Molecule 

CHF 

SDQ-MBPT(4) 

Experiment 

HF 

p 

0.758 

0.709 

0.707 

a 

4.89 

5.58 

5.52 

h2o 

p 

0.784 

0.735 

0.724 

a 

8.53 

9.54 

9.64 

nh3 

P 

0.614 

0.578 

0.579 

s 

1  .1 

14.2 

14.8 

Li  F 

u 

2.56 

2.49 

2.49 

ct 

7.48 

9.83 

•  •  • 

BeO 

p 

2.96 

2.52 

•  •  • 

T5T 

20.9 

29.73 

. . . 
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experiment  is  poor.  In  terms  of  xj^>  our  best  CHF  calculation  predicts 
-70.8,  while  for  basis  A,  we  obtain  -85.8.  These  are  improvements  over  other 
CHF  calculations,  but  obviously,  even  with  superior  basis  sets,  there  are 
large  effects  that  contribute  to  the  hyperpolarizabilities  that  are  not  yet 
included. 

The  largest  single  effect,  we  believe,  is  the  effect  of  electron 
correlation.  This  is  demonstrated  in  Table  II.  In  our  best  basis  set,  B, 
correlation  causes  a  factor  of  four  increase  in  szzz  with  a  more  modest 
increase  in  BZyy-  Tbe  increase  in  basis  A  is  not  as  great,  but  still  amounts 
to  a  factor  of  more  than  two  in  Qzzz-  It  is  not  entirely  clear  why  such  large 
effects  are  observed  due  to  correlation,  but  it  apparently  pertains  to  a  much 
more  accurate  description  of  the  long-range  part  of  the  charge  density  that 
is  being  preferentially  sampled  by  properties  like  hyperpolarizabilities. 

Also,  the  lone  pair  of  electrons  on  the  z-axis  is  expected  to  be  quite 
sensitive  to  correlation  effects. 

In  Table  III  a  comparison  of  NHj  with  HF  and  H2O  is  presented. 
Correlation  changes  8  by  almost  100  percent  of  the  CHF  result  for  NH-j,  while 
being  respectively  28  percent  and  50  percent  for  HF  and  H2O.  This  may  be 
due  to  the  single  lone  pair  in  NH3  contributing  its  entire  effect  in  one 
direction  while  the  two  pairs  in  HgO  and  three  in  HF  have  somewhat  opposing 
effects.  Whatever  the  reason,  with  correlation  our  prediction  for  xj^  for 
NHg  is  -138.7  XI 0~ 3 3  esu/molecule  in  much  better  agreement  with  the  dc- induces 
second  harmonic  prediction  of  -209  +5xl0"33  esu/molecule  than  all  previous 
calculations  as  shown  In  Table  III.  It  is  apparent  that  the  agreement  can 
still  be  Improved,  but  it  is  certainly  evident  that  only  a  correlated  approach 
has  any  hope  of  successfully  predicting  hyperpolarizabilities. 


TABLE  II 
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a  Experimental  results  are  those  of  Miller  and  Ward,  Reference  14. 
b  D.  P.  Sa  try,  Chem.  Phys.  Letters  53,  568  (1978),  CHF  calculations. 
c  S.  P.  Liebman  and  J.  W.  Moskowitz,  J.  Chem.  Phys.  54,  3622  (1971),  UCHF  calculations. 
d  Reference  ,  CHF  calculations. 
e  Reference  ,  CHF  calculations. 
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The  components  of  the  y  hyperoolarizability  are  similarly  observed 
to  be  greatly  affected  by  basis  set  and  by  correlation.  The  behavior  is 
rather  more  predictable  than  in  the  case  of  however,  since  improving  the 
basis  set  increases  the  size  of  and  correlation  further  increases  it. 

The  correlation  effect  on  y22Z2  is  72  percent  of  the  CHF  result,  but  even 
the  more  diffuse  basis  contributes  this  large  a  percentage  change.  3  on 
the  other  hand,  tends  to  become  smaller  with  improvements  in  basis  sets, 
while  becoming  larger  due  to  correlation  effects.  Since  the  correlation 
effect  is  usually  larger  in  a  better  basis,  though,  the  basis  set  effect 
is  partially  offset  when  fully  correlated  calculations  are  made. 

The  number  in  parenthesis  in  Table  1  show  the  change  incurred  when 
a  smaller  field  strength  of  0.005  instead  of  0.01  is  used.  This,  indicates 
the  degree  of  contamination  remaining  in  the  finite  field  method.  Although 
there  is  almost  a  one  unit  change  in  3222>  it  is  clear  that  contamination  will 
have  little  net  effect  on  the  observed  results. 
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V.  REVIEW  OF  RESEARCH  ACCOMPLISHMENTS 


Linder  Air  Force  Office  of  Scientific  Research  sponsorship,  a  number 
of  accomplishments  have  been  made  in  this  study  that  we  believe  make  a  signifi¬ 
cant  contribution  to  the  problem  of  atomic  and  molecular  hyperpolarizabilities, 
and  the  related  nonlinear  optical  effects  of  gases. 

A.  This  project  provides  the  first  studies  of  molecular 
hyperpolarizabilities  ever  attempted  which  include 
electron  correlation.  Molecules  studied  include,  N^, 

HF,  H20,  NH3,  Bg0 j  and  Li F. 

B.  This  work  has  demonstrated  that  correlation  effects 
are  tremendously  important  in  a  predictive  theory  of 
hyperpolarizabilities.  In  the  case  of  NHj,  Be0  and  Li F , 
correlation  can  change  the  hyperpolarizability  by  more  than 
a  factor  of  two. 

C.  To  make  state-of-the-art  quantum  mechanical  calculations 

of  hyperpolarizabilities,  the  theory  of  many-body  perturbation 
theory  (MBPT)  has  been  extended  to  permit  the  inclusion  of 
all  important  correlation  effects  in  hyperpolarizabil ity 
calculations. 

D.  A  group  of  highly  efficient  ab  initio  computer  programs  based 
upon  MBPT,  and  that  can  use  quite  large  basis  sets,  have 
been  developed  to  make  possible  reliable  calculations  of 
hyperpol ari zabi 1 1 ti es . 

-  ; 

_  t 
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E.  We  have  also  generalized  the  calculation  to  include 
frequency-dependence  and  to  treat  the  important  correlation 
effects  in  an  infinite-order  correlated  treatment  based 
upon  the  coupled-cluster  formulation. 

F.  For  the  first  time,  reasonable  agreement  between  theoretical 
calculations  and  second-harmonic  generation  experiments 

for  6=  (6zz2+8xx2+SyyZ )  has  been  achieved.  For  H20,  Ward 
and  Miller  report  yi*)  to  be  -  94+4  x  1C  3'esu/molecule, 
while  our  predictions  give  -72  and  -80  x  1033  esu/molecule 
for  the  static  hyperpolarizability. For  NH3  experiment  reports 
-209  x  10  33esu/molecule  while,  our  calculations  give  at 
w=0,  138.7  x  IQ'33  esu/molecule.  All  previous  studies  of  NH3  were 
in  much  poorer  agreement  with  experiment,  often  even  getting 
an  incorrect  sign. 

G.  In  another  investigation,  which  is  also  the  first  of  its  kind, 
we  have  studied  the  B  hyperpolarizability  tensor  as  a  function 
of  geometry  for  H,0.  It  is  found  that  6  is  extremely  sensitive 

C  *\, 

to  molecular  vibrations,  with  its  value  changing  by  nearly 
50  percent  just  from  equilibrium  to  the  maximum  extent  of  the 
zero  point  vibration,  for  the  symmetric  stretch  mode  (See 
Appendix  B). 

H.  From  (G),  it  is  apparent  that  proper  account  of  the  sensi¬ 
tivity  of  hyperpolarizabilities  with  vibration  is  important 
if  theoretical  and  experimental  numbers  are  to  coincide. 


Additional  studies  have  been  made  to  determine  the  types  of 
basis  sets  required  to  obtain  reliable  higher  polarizabilities. 
This  is  extremely  critical  to  the  theoretical  predictions. 

A  comparison  with  completely  numerical  coupled  Hartree-Fock 
studies  of  the  HF  molecule  were  undertaken,  since  these 
results  should  entail  essentially  no  basis  set  error.  The 
types  of  basis  sets  used  are  found  to  agree  with  the  fully 
numerical  CHF  calculations  to  within  'vIO  percent  (See  Appendix  A). 
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Molecular  hyperpolarizabilities.  I.  Theoretical  calculations  including  correlation 

Rodney  J.  Bartlett  and  George  D.  Purvis  III 
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Static  polarizabilities  and  hyperpolarizabilities  for  molecules  are  investigated  at  the  correlated  level.  The 
(inite-field,  coupted  Hartree-Fock  theory  is  used  as  a  zeroth-order  approximation,  with  correlation  included 
by  using  the  linked-diagram  expansion  and  many-body  perturbation  theory,  that  includes  single,  double,  and 
quadruple  excitation  diagrams.  The  theory  is  illustrated  by  studying  the  hydrogen  fluoride  molecule.  It  is 
demonstrated  that  the  correlation  effect  for  the  hyperpolarizabilities  if  and  ?  can  be  quite  large.  The 
average  polarizability  and  dipole  moment  of  HF  are  in  excellent  agreement  with  experiment.  The  relative 
importance  of  the  various  types  of  diagrams  contributing  to  electric  field  properties  are  discussed.  The 
dependence  of  the  computed  hyperpolarizability  on  basis  sets  is  also  investigated. 


1.  INTRODUCTION 

The  nonlinear  optical  properties  of  gases  have 
been  of  experimental  interest  for  several  years.1"7 
Recently,  a  number  of  real  and  potential  laser 
applications  utilizing  the  frequency-tripling  pro¬ 
perties  of  metal  vapors  and  other  atomic  gases  have 
been  suggested.8'10  Similar  devices  employing 
polar  molecules  oriented  in  an  electric  field  can 
be  envisioned  for  frequency-doubling  applications. 
The  utility  of  such  novel  devices  will  be  ultimately 
determined  by  properties  of  the  atomic  and  mole¬ 
cular  higher  polarizability  tensors,  but  little  in¬ 
formation  on  these  quantities  currently  exists. 

Hyperpolarizabilities  may  be  deduced  experi¬ 
mentally  from  direct  observations  of  harmonic 
generation  in  gases3"5  and  from  measurements  of 
the  Kerr  effect.8’7  However,  the  experiments  are 
difficult  and  the  range  of  uncertainty  is  often  large. 
Since  hyperpolarizability  is  a  property  of  a  single 
atom  or  molecule,  it  may  also  be  predicted  from 
quantum- mechanical  calculations.  Hence  it  is 
important  to  have  reliable  theoretical  predictions 
for  higher  polarizabilities  to  complement  the 
experimental  efforts. 

While  some  coupled  Hartree-Fock  (CHF)“  level 
calculations  of  atomic  hyperpolarizabilities  are 
available,  mostly  for  inert  gases, 13  very  few  ab 
initio  calculations  of  molecular  hyperpolarizabili¬ 
ties  have  been  attempted,  ‘I"1*  and  none  of  these 
has  yet  shown  any  kind  of  agreement  with  experi¬ 
ment8  The  molecular  calculations  which  have 
been  made  frequently  employ  inadequate  basis  sets 
or  are  carried  out  at  the  level  of  uncoupled  Har¬ 
tree-Fock  perturbation  theory.13*17  Only  a  few 
previous  studies  have  even  used  the  full  CHF 
method,  18"‘*  with  no  work,  at  all,  on  molecular 
hyperpolarizabilities  at  the  correlated  level.  Re¬ 
cent  papers  have  addressed  the  question  of  picking 
adequate  basis  sets  for  polarizabilities,  1*"ao  but 
the  accurate  Inclusion  of  correlation  effects  is 
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still  expected  to  be  a  necessity  if  a  reliable  pre¬ 
dictive  theory  of  hyperpolarizabilities  is  to  be 
developed.  Some  recent  communications  have 
demonstrated  the  importance  of  correlation  for 
the  dipole  polarizability,18’20"22  and  one  would  anti¬ 
cipate  an  equally  large,  or  even  larger,  correla¬ 
tion  effect  for  the  hyperpolarizabilities. 

Recent  developments  in  many-body  perturbation 
theory  (MBPT)23*27  and  the  coupled-cluster  ap¬ 
proaches  (CCA),  25’28’29  have  made  it  possible  to 
include  correlation  in  a  sufficiently  tractable 
manner  to  include  a  very  large  part  of  the  net 
correlation  effect  in  molecular  calculations25’30 
even  though  comparatively  large  basis  sets  are 
required  in  hyperpolarizability  determinations. 

In  the  following,  we  report  a  study  of  the  hyper¬ 
polarizabilities  of  the  HF  molecule,  considering 
correlation  effects  due  to  single,  double,  and 
quadruple  excitations.25  HF  provides  a  convenient 
example  for  this  initial  study  since  a  numerical 
CHF  result  has  been  obtained  for  the  parallel  com¬ 
ponents  of  the  second-  and  third-order  polariza¬ 
bility  by  Christiansen  and  McCullough, 18,19  and 
this  result  may  be  used  to  eliminate  some  uncer¬ 
tainty  in  chosing  a  reliable  basis  set. 

II.  HIGHER  DIPOLE  POLARIZABILITIES 

The  energy  of  a  molecule  in  an  external  field 
$  may  be  written  in  a  power  series  as 

W(J)  =  W<0>-  H0, -WY' 

X  A  tfi t"  ’  '  ■ 

where  the  summation  over  repeated  indices  is 
assumed.  6  f  is  a  component  of  the  external  field. 
lit  is  a  component  of  the  permanent  dipole,  au  a 
polarizability,  and  and  yllt,  are,  respectively, 
the  second  and  third  polarizabilities.'"8  A  com¬ 
ponent  of  the  total  dipole  moment,  />,,  is  obtained 
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from  the  derivative  dW  'dS ,,  which  is  composed 
of  the  permanent  moment  and  an  induced  moment 
due  to  the  higher-order  terms. 

(iW  \ 

Jg-J  =  -Pi  =  A  *  (2!  A 

+  (V.)-lYinlSjS  >&,-■■  ■  (2) 

The  production  of  second  and  third  harmonic 
generation  follows  from  the  fact  that  the  7 and  7 
in  term  in  Eq.  (2)  involve  products  of  electric 
field  components.  Hence  for  S ,  =S0I  siniuf,  we 
have 

S 2,  =  <Sj0  sin2cuf  =  |<S  |0(1  -  cos2uf) , 

and  the  cos2u)f  term  gives  rise  to  frequency  doub¬ 
ling.1  Similarly,  the  third  power  in  S  is  responsi¬ 
ble  for  frequency  tripling.1 

The  polarizabilities  o’,  j3,  T are  said  to  be  se¬ 
cond-,  third-,  and  fourth-order  properties, 
respectively,  based  on  the  order  of  the  external 
field  in  the  energy  expansion  of  Eq.  (1),  while  the 
permanent  dipole  moment  is  first  order.  It  may 
be  shown  that  a  Hartree-Fock  wave  function  will 
predict  a  first-order  property,  like  the  per¬ 
manent  dipole  moment,  to  a  comparatively  good 
accuracy  because  the  first  correlation  corrections 
vanish, 31,32  but  there  is  no  similar  reason  to  ex¬ 
pect  the  CHF  theory11*33  to  be  adequate  lor  the 
higher-order  properties.1** 21,22  Furthermore,  for 
even  highly  accurate  first-order  properties  and 
especially  higher-order  properties,  correlation 
must  be  considered. 

in.  MANY-BODY  PERTURBATION  THEORY 

The  coefficients  in  Eq.  (1)  can  be  obtained  direct¬ 
ly  from  perturbation  theory  or  by  using  finite- 
field  methods.  Choosing  the  latter,  temporarily, 
we  may  write  the  Hamiltonian  for  the  perturbed 
molecule  as 

x(S)=x+\a(6)=F0  +  v  +  \n(S)  (3) 

and 

F(3)=F0  +  \n(S),  (4) 

xn =  ou(i)  =  - (5) 

rW-E  r-J-E  *«,?).  (8) 

«</  i 

F0  is  the  usual  independent  particle  self-consistent- 
field  (SCF)  effective  Hamiltonian;  hence 

F(S)’^lh(i)+u(i,S).S-r(i)]  (7) 

i 

u(  1,J)-E  f  (8) 

i  J 


and  in  the  general  case  the  molecular  orbitals 
{*,}  are  field  dependent. 

If  we  are  only  interested  in  the  SCF  solution  in 
the  presence  of  the  field,  then  V  in  Eq.  13)  may  be 
neglected.  This  defines  the  CHF  model,  with  the 


field-dependent  solutions 

F(J)C,(J)  =  A£,l£  )<,{£).  (9) 

The  molecular  orbitals  \P(S)  are  defined  as  an 
expansion  in  terms  of  an  atomic  orbital  basis  set 

!«>: 

X,($)  =  !!>£,(?).  (10) 

From  Eqs.  (6)-(9)  it  follows  that 

F(S)*0(S)=  W0(S)<P0(S),  (11) 

where 

<i.a(<f)=  aU,(D-  •  •  x„M.  (i2) 

) «  W0&)  *  <*0  |  V$)  |*0> .  (13 ) 


In  the  CHF  case,  the  SCF  equations  are  solved 
for  fixed  (finite  field)  values  of  the  field  strength 
S.  From  these  solutions,  the  CHF  energy  may  be 
obtained  as  a  function  of  S  as 


-§£<x,(S)x,(S)llx,(£>x,(?»-  (14) 

C7-i 

The  double-bar  notation  designates  the  two-elec¬ 
tron  integral, 

<X»X,IIXrX,> 

=  /  dr,  /  dT2X;(l)x,*(2)r-1(l  -PJ 

*X,U>X,(2).  (15) 

The  CHF  polarizabilities  and  hyperpolarizabilities 
may  be  derived  from  second  and  higherderivatives 
of  Wclir(3)  with  respect  to  the  field  strength. 
Alternatively,  the  induced  dipole  moment  can  be 
calculated  from  the  wave  function,  Eq.  (12).  with 
the  polarizability  and  hyperpolarizabilities  ob¬ 
tained  as  first,  second,  and  higher  derivatives 
of  the  induced  dipole  moment.  The  rigorous 
equivalence  between  the  dipole  procedure  and  the 
energy-based  procedure  depends  upon  the  satis¬ 
faction  of  the  Hellman-Feynman  theorem  for  the 
approximate  wave  functions.22 

As  an  alternative  to  the  CHF  model,  a  perturba¬ 
tion  expansion  of  F{§)  and  u(S )  in  powers  of  J  may 
be  made.  This  procedure,  subject  to  orthonorma- 
Uty  of  the  molecular  orbitals,  leads  to  the  coupled 
perturbed  Hartree-Fock  (CPHF)  method.34'3*  The 
results  of  CPHF  and  CHF  are  formally  the  same.11 
In  order  to  go  beyond  the  CHF  level  and  include 
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FIG.  1.  Single-particle  corrections  to  the  energy  dia¬ 
grams.  ( — *)  is  the  negative  of  the  Fock  potential.  If 
occupied  and  excited  one-particle  orbitals  are  deter¬ 
mined  as  SCF  solutions  at  any  field  strength,  then  these 
diagrams  component  mutually  cancel. 


effects  of  correlation,  the  perturbation  V{S)  of 
Eq.  (6)  needs  to  be  considered  in  more  detail.  In¬ 
cluding  this  perturbation,  the  energy  for  the  per¬ 
turbed  molecule  can  be  obtained  from  the  linked- 
diagram  expansion24  as 

w(S)  =  wCHT(S) 

+£  <<f>o|A5)[W'0(?)-F(?)]-‘V(J)]»|<f>a>f 


The  terms  defined  by  the  summation  in  Eq.  (16) 
are  said  to  be  the  correlation  corrections  to  the 
energy  for  a  given  field  strength.  Just  as  in  the 
zero-field  case,  by  solving  the  SCF  equations  in 
the  presence  of  the  field  as  in  Eqs.  (8)— (10),  the 
SCF  cancellations  shown  in  Goldstone-diagram 
form  in  Fig.  1  ( —  x  is  the  negative  of  the  Fock 


£  c 

IA)  iB>  .C 

<2^  dr  0*r 


FIG.  2.  Zeroth-  and  first-order  relaxation  diagrams 
for  a  second-order  property.  These  terms  are  all  in¬ 
cluded  in  a  CHF  calculation  for  a  second-order  propertv 


potential)  are  still  maintained,  so  that  no  diagrams 
containing  these  parts  are  required. 

The  coefficients  in  Eq.  (1)  are  given  as  deriva¬ 
tives  of  K'(<? )  with  respect  to  field  strength.  In 
the  finite-field  approach,  U '(3)  must  be  computed 
for  several  field  strengths  to  allow  the  determina¬ 
tion  of  the  polarizabilities,  which  requires  taking 
small  differences  between  large  terms.  If  the 
differentiation  is  made  initially,  separate  hnked- 
diagram  expansions  for  p,  or,  d,  >.  etc.,  in  terms 
of  the  zero-field  SCF  orbitals  [x/01-xi-  «,(0) 

=  «^],  can  be  used  to  evaluate  the  small  polariza¬ 
bilities  directly. 

In  the  latter  case,  ordinary  double  perturbation 
theory  gives 

*•0 

xK»'i,-^)',^(V  +  An-AlF')]»|+a)  (17) 

and.  for  some  component  St, 

W(«,).l»r,.(Wu0*Hr|tI.U'1.1.-  ■  •  >5, 

U8) 

=  lV0.-  *  «■•••• 

(19) 

Considering  a  second-order  property  for  illus¬ 
tration,  with  the  definition 

*„=  PWa-Fn)-'P,  (20) 

where  P  is  the  projector  for  the  orthogonal  com¬ 
plement  to 


(A,.  (£,)  *  K,  =  (*0  l+o)  ♦  (*o  -  Wl%0)RaV  t*  J  *  <*, 1  Vfi0(x  -  lFu0Mlo*  l*o> 

+  <*ok«0(V-lF0fl)flffx|*0>+-  •  •  . 


The  first  term  on  right  in  Eq.  (21),  is  given 

by  the  uncoupled  Hartree-Fock  (UCHF)11  approxi¬ 
mation  [Fig.  2(A)],  which  is 


,  T  T 

)  L*  L*t  • 

((oee)<(u«e)  M 


while  the  higher  terms  involve  the  field-indepen¬ 
dent  perturbation  F(0).  The  Goldstone  diagrams 


for  second  order  in  an  external  field  and  for  zeroth 
and  first  order  in  the  two-electron  perturbation 
F(0)  are  shown  in  Fig.  2.  The  symbol  ( — •)  indi¬ 
cates  interaction  with  the  external  field. 

It  is  important  to  recognize  that  even  though 
diagrams  (B)-(G)  of  Fig.  2  involve  a  two-electron 
vertex,  these  diagrams  and  selection  oi  others 
in  all  higher  orders  in  i'(0)  are  included ^T’,4  in  the 
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CHF  results  or  in  the  perturbation-theory-equiva- 
lent  CPHF.'%'“  At  least  In  the  static  case  these 
terms  should  not  be  considered  as  correlation 
corrections.  These  terms  arise,  instead,  from 
the  effective  Fock  potential  uti ),  now  written  in 
terms  of  perturbed  orbitals,  which  reflects  the 
relaxation  of  the  orbitals  to  the  external  perturba¬ 
tion.'1  Beyond  first  order  in  l\0).  additional 
diagrams  arise  which  are  not  included  in 
CPHF.  M’ ,T  and  these  terms  constitute  true  corre¬ 
lation  corrections.  In  the  finite-field  approach, 
all  terms  which  arise  from  derivatives  of  the 
summation  in  Eq.  (16)  are  actual  correlation  cor¬ 
rections.  while  the  initial  CHF  calculation  pro¬ 
vides  the  zeroth-  and  first-order  energy  correc¬ 
tions,  subject  to  the  field-dependent  orbitals. 
(CPHF  is  the  static  equivalent  of  time-dependent 
Hartree-Fock  theory  or  the  random-phase  approxi¬ 
mation.  which  is  similarly  recognized  to  sum 
selections  of  many-body  propagator  diagrams  to 
all  orders. 38,39  In  the  time-dependent  context,  the 
terms  so  summed  are  frequently  referred  to  as 
“dynamic"  correlation.39) 

Even  though  the  direct  determination  of  the 
polarizabilities  is  a  distinct  advantage,  it  offers 
the  disadvantage  that  separate  sets  of  programs 
must  be  written  to  evaluate  the  diagrams  for  each 
order  in  the  external  perturbation,  while  the 
finite-field  method  allows  all  polarizabilities  to 
be  obtained  from  simply  executing  “enough”  energy 
calculations  at  a  series  of  field  strengths.  This 
also  permits  one  to  exploit  the  theories  and  pro¬ 
grams  that  have  been  developed  for  the  usual 
correlation  problem25  which  typically  offer  a  much 
more  sophisticated  level  of  treatment  for  the 
correlation  than  would  be  convenient  to  develop 
for  each  individual  order  in  an  external  perturba¬ 
tion.  Also,  the  dichotomy  into  correlation  and 


relaxation  effects  is  transparent  :  >r  the  finite-field 
method. 

IV  M  MERIC \L  RESI  LTS 

In  the  finite-field  methods,  it  is  necessary  t  ' 
obtain  the  various  polarizability  components  from 
formulas  for  the  energy  or  dipole  moment.  In  the 
present  work.  either_procedure  may  be  used  at  the 
CHF  level,  since  <t> , t<S  )  satisfies  the  Hellman- 
Feynman  theorem,  but  derivatives  of  IT-,6’  i  will  be 
used  for  the  correlation  corrections.  By  consider¬ 
ing  Eq.  (1).  it  may  be  shown  that  the  finite-diffe¬ 
rence  formulas  listed  in  T able  I  hold  for  the  various 
polarizabilities.  Each  of  these  formulas  is  obtain¬ 
ed  by  excluding  all  even  or  odd  terms  in  Eq.  (1)  by 
using  positive  and  negative  field  strengths  of  the 
same  magnitude.  This  ensures  that  the  contamina¬ 
tion  from  the  next  higher  term  in  the  power  series 
is  completely  removed,  leaving  only  the  next  high¬ 
er  term  of  the  same  type  (i.e.,  even  or  odd),  which 
is  about  four  orders  of  magnitude  smaller  (at  the 
field  strengths  used  here)  than  any  polarizability 
being  determined.  Hence,  essentially  no  error  is 
anticipated  from  higher  polarizability  contamina¬ 
tion.  At  the  same  time,  however,  field  strengths 
of  adequate  size  must  be  used  to  ensure  that 
significant  energy  differences  are  obtained.  In 
this  work  field  strengths  of  0.0,  0.01.  and  0.02  a.u. 
are  found  to  be  suitable. 

The  formulas  in  Table  I  are  general,  but  if  some 
symmetry  is  present,  formulas  1.3  and  1.6  be¬ 
come  much  simpler.  If  the  molecule  has  a  rotation 
axis,  with  i  representing  the  direction  of  this  axis, 
and;  is  perpendicular  to  i,  group-theoretical  con¬ 
siderations  show  that  the  components.  p,„,  o(/, 
yUjl,  and  ylllt  are  vanishing  at  zero  field  strength. 
This  results  in  particularly  simple  formulas  for 


TABLE  I.  Energy  formulas  for  finite-field  calculations  of  dipole  moments,  polarizabilities,  and  hvperpolariz- 
abllities. 1 


Odd  order 

utS  »(-Si)|  +  ^l»'<2<S,)_»'<-2Sl>j+0<6>  (11 

Even  order 

all$i=+$ w<0>  — f  IW(«S,)+  »(_$,))  +  fo[WU&,)*  W(_2(S,»+0(D  <41 

r„,j(Sf=+4{R'((S,)+  »■(-<£  ())-[  W  (2<S,)+  R'(-25j)|-6H'(0)+O(€)  (3) 

+  [WS;)*  w{-S,)]-2wm+o<(  i 


For  a  molecule  with  a  rotation  axis  in  the  direction  «,  atj,  y jtu,  and  ytllj  will  be  zero  by  symmetry. 
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FIG.  3.  All  second-  and  third-order  correlation  dia¬ 
grams  provided  that  SCF  (CHF)  orbitals  are  used  for 
the  occupied  and  excited  one-particle  states.  Anti¬ 
symmetrized  vertices  are  assumed  in  these  diagrams. 

the  nonvanlshing  components. 

The  finite-field  correlation  corrections  included 
are  shown  as  antisymmetrized  diagrams  in  Figs. 
3-6,  where  it  is  understood  that  the  hole  and 
particle  lines  are  field  dependent  (or  "dressed”) 
and  that  finite  differences  of  these  quantities  must 
be  taken  before  the  polarizabilities  are  obtained. 
(See  Refs.  27  and  40  for  the  rules  and  algebraic 
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FIG.  4  All  single-excitation  fourth-order  diagrams 
provided  that  SCF  (CHFi  orbitals  are  used  for  the 
occupied  and  excited  one-particle  states.  Antisymme¬ 
trized  vertices  are  assumed. 

formulas  corresponding  to  these  antisymmetrized 
diagrams. ) 

These  diagrams  consist  of  all  terms  that  occur 
through  fourth-order  in  the  correlation,  subject 
to  field-dependent  SCF  orbitals  that  would  arise 
due  to  configuration-interaction  (Cl)  single,  double, 
and  quadruple  excitations.35  Triple  excitations 
also  contribute  in  the  fourth-order  energy,  but 
these  are  excluded. 

It  is  well  known  that  Cl-type  single  excitations 
are  quite  important  in  determining  properties  other 
than  the  energy,  since  operators  such  as  &  •  r.  will 
mix  single  excitations  directly  with  an  SCF  unper¬ 
turbed  wave  function.  In  the  present  work,  the 
predominant  effect  of  single  excitations  is  intro¬ 
duced  by  means  of  the  initial  CHF  calculations, 
which,  as  described  in  Sec.  III.  is  responsible  for 
summing  a  series  of  diagrams  involving  such 
single-excitation  vertices  as  occur  in  Fig.  2  (i.e. , 

A T*)  to  all  orders.  This  “dresses"  the  hole  and 
particle  lines  involved  in  Figs.  3-6.  The  remain¬ 
ing  single  excitations  accounted  for  by  Fig.  4  in¬ 
volve  the  smaller  effect  of  single  excitations  of 
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FIG.  5.  All  double-ex¬ 
citation  fourth-order  dia¬ 
grams  provided  that  SCF 
(CHF)  orbitals  are  used 
for  the  occupied  and  ex¬ 
cited  one-particle  states. 
Antisymmetrized  vertices 
are  assumed. 
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FIG.  6.  All  quadruple-excitation  fourth-order  dia¬ 
grams,  provided  that  SCF  (CHF)  orbitals  are  used  for 
occupied  and  excited  one-particle  states.  Unlinked  dia¬ 
grams  (H)  and  (I)  do  not  contribute  to  the  correlation 
energy.  Antisymmetrized  vertices  are  assumed. 


dressed  lines  that  enter  into  the  fourth-order 
energy  only  by  interacting  through  dressed  double 
excitations.  The  very  high-order  coupling  of  the 
first  category  of  single-excitation  effects  with 
correlation  due  to  double  and  quadruple  diagrams 
by  using  "dressed”  hole  and  particle  lines  is  one 
of  the  additional  advantages  of  the  finite-field- 
based  methods. 

The  basis  set  for  the  present  study  of  HF  is 
chosen  essentially  following  the  prescription  of 
Christiansen  and  McCullough  (CM).19  We  start 
from  Dunning’s11  5s3/>  contraction  of  Huzinaga's12 
9s5/>  basis  set  for  flourine,  and  Dunning’s  3s  con¬ 
traction  for  hydrogen  Slater  exponent  1.2.  These 
functions  are  augmented  by  a  two-Gaussian  fit  to 
a  d-STO  for  F  with  exponent  3.358,  and  a  two- 
Gaussian  fit  to  ap-STO,  exponent  2.082  for  H.13 
This  results  in  a  5s3plrf/3sl/>  initial  basis  set. 
This  basis  is  then  augmented  in  a  "well-tempered" 


.’(> 

fashion.  The  smallest  exponent  for  each  l  quantum 
number  per  atom  is  multiplied  by  0.4  to  obtain  a 
new,  diffuse  function  to  assist  in  describing  the 
polarizability.  The  resulting  basis  6<4/>2d  4s2/> 
is  expected  to  be  adequate  to  describe  the  dipole 
polarizability.  To  also  attempt  to  account  f  jr  the 
higher  polarizabilities,  the  new  exponent  is  again 
multiplied  by  0.4  to  add  another  />  and  d  function  on 
F.  and  an  s  and  p  function  on  H.  Finally,  the  most 
diffuse  d  exponent  on  F  was  again  multiplied  by 
0.4  and  this  function  added  to  the  basis  set.  We 
also  consider  the  addition  of  an  extra  diffuse  s 
function  on  F  and  the  exclusion  of  the  most  diffuse 
s  and  p  on  H  and  the  d  function  on  fluorine.  This 
provides  four  different  basis  sets. 

A  GTO  basis  set  similar  to  these  has  been  shown 
to  account  for  the  numerical  CHF  results  for 
oi ^ ,  d,„.  and  to  within  a  few  percent.'1  How¬ 
ever.  those  authors  were  unable  to  reduce  the  er¬ 
ror  in  the  CGTO  value  for  below  7^.  When  we 
attempted  to  use  their  basis  set.  with  all  six  possi¬ 
ble  Cartesian  components  of  the  d  functions,  xx. 
yv ,  zz,  xz.yz,  xy  (CM  use  only  the  five  normal  d- 
function  components),  we  found  an  approximate 
linear  dependency  which  made  our  correlated  cal¬ 
culations  unstable.  We  were  able  to  resolve  this 
problem  and  also  substantially  reduce  the  error 
in  the  SCF  ji,gt  component  by  choosing  a  tighter 
rf-polarization  function  on  F  (STO  t,  =  3.358  com¬ 
pared  to  t,  =  2.25  for  CM)  and  then  adding  an  addi¬ 
tional  diffuse  d  function  to  the  F  atom.  Both  the 
tight  and  diffuse  d  functions  seem  to  be  important 
as  illustrated  in  Table  II.  Using  a  tight  d  function 
while  excluding  the  most  diffuse  rf-function  (bases 
A  and  B)  results  in  a  dtu  that  differs  from  the 
numerical  result  by  ll^.  In  the  CM  basis,  which 
has  a  very  diffuse  d  function  but  excludes  a  tight  po¬ 
larization  d  function,  a  similar  error  in  d,„  occurs. 
Bycombiningboth.  as  inbasesC  andD,  thiserroris 
greatly  reduced. 

The  difference  between  bases  C  and  D  is  the 


TABLE  n.  Comparison  of  CHF  results  for  HF  in  various  basis  sets  with  numerical  CHF. 
(R»  1.7328  bohrs;  values  In  a.u.) 


Basis  set  (GTO) 

WTO) 

u* 

y**u 

CM  <6s5/>3d/4s3p)6 

-100.0535 

0.759 

5.80 

-9.0 

310 

A  (6s5p3d/5s3p)z 

-100.0563 

0.757 

5.72 

-9.3 

250 

B  as5p3d/5s3p)c 

-100.0563 

0.757 

5.72 

-9.2 

260 

C  (6s5p4d/4s2p)c 

-100.0563 

0.759 

5.73 

-8.4 

250 

D  (6s5p4d/5s3p)c 

-100.0565 

0.758 

5.76 

-8.5 

2S0 

Numerical  CHF6 

-100.0706 

0.756 

5.76 

-8.3 

320 

*  Coordinate  system  is  chosen  such  that  dipole  moment  (F‘H#)  is  positive. 

4  Results  of  Ref.  16. 

c  Basis  sets  A,  B,  C,  and  D  use  all  six  Cartesian  Gaussian  d -orbital  functions  xx,  vy,  «, 
xy,  vs,  and  xx. 
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TABLE  IV.  Contributions  of  different  diagrams  to  dipole  moment,  polarizability,  and  hyperpolarizabilities  of  HF.  (Basis  set  U;  K=  1.7328  bohrs;  values  in 
atomic  units.*) 
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exclusion  of  the  most  diffuse  $  and  p  functions  on 
H.  This  is  seen  to  have  almost  no  effect  on  p  while 
having  a  somewhat  larger  effect  on  >„„■  The 
observation  that  the  diffuse  H  s  and  p  functions 
have  only  a  small  effect  is  partially  due  to  the  fact 
that  the  diffuse  region  of  the  charge  cloud  is  al¬ 
ready  accounted  for  by  the  other  highly  diffuse 
functions  on  F.  However.  y„u  is  usually  larger 
the  more  diffuse  functions  that  are  included  in  the 
basis,  as  supported  by  the  results  of  bases  sets 
C  and  D. 

It  is  clear  from  Table  II  that  good  agreement 
with  the  numerical  results  for  p.,.  a„.  P,„.  and 
y„„  can  be  achieved  with  the  present  basis  sets. 

It  does  not  necessarily  follow  that  the  other  com¬ 
ponents  of  3,  2.  and  y  are  as  well  described  or 
that  this  basis  is  entirely  adequate  for  the  cor¬ 
relation  corrections  to  the  polarizabilities,  but  the 
good  agreement  with  the  numerical  results  for 
the  parallel  components  at  least  provides  an  indi¬ 
cation  that  the  final  correlation  corrections  to  the 
different  polarizabilities  should  be  indicative  of 
the  size  of  the  true  correlation  corrections,  which 
is  the  primary  objective  of  this  study. 

In  the  present  correlated  calculations,  all  the 
single,  double,  and  quadruple  excitation  diagrams 
that  arise  through  fourth  order  in  the  l'(<?)  per¬ 
turbation  are  included.  This  model  willbe  referred 
toasSDQ-MBPT(4).  These  diagrams  are  evaluated 
at  various  field  strengths,  from  which  the  equa¬ 
tions  of  Table  I  are  used  to  provide  the  different 
polarizabilities. 

To  carry  out  such  finite-field  calculations  suc¬ 
cessfully,  it  is  necessary  to  ensure  about  eight- 
decimal-place  accuracy  in  all  computations,  from 
the  initial  molecular  integrals  to  the  actual  dia¬ 
gram  evaluation.  The  current  computations  em¬ 
ploy  the  MBPT  program  system  developed  at 
Battelle.44  Results  for  the  effect  of  correlations 
on  various  properties  are  shown  in  Table  III. 

The  6s5p4d/5s3p  basis  set  is  seen  to  be  capable 
of  providing  about  78%  of  the  observed  field-free 
valence-shell  correlation  energy.  This  is  con¬ 
sistent  with  other  calculations,  where  we  have 
shown  that  a  5s3pld /3slp  basis  typically  accounts 
for  three-quarters  of  the  valence-shell  correlation 
energy.25'45  since  the  remaining  functions  in  the 
current  bases  are  generally  too  diffuse  to  contri¬ 
bute  much  to  the  field-free  correlation  effect  On 
the  other  hand,  correlation  involving  the  normal 
and  diffuse  functions  is  important  in  the  presence 
of  the  field,  so  some  balance  between  a  reasonable 
description  of  the  valence-shell  correlation  and 
the  long-range  tails  of  orbitals  is  still  expected  to 
be  significant  in  obtaining  good  correlated  results 
for  polarizabilities. 

The  effects  of  correlation  on  the  other  properties 
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listed  in  Table  11!  vary  from  about  7  7  for  the  dipole 
moment  to  12"  for  the  polarizability,  and  a  rather 
dramatic  change  of  about  22"  for  the  d  hyperpolar¬ 
izability.  A  similarly  large  change  is  observed 
for  the  y„„  component  of  the  hyperpolarizability 
of  307-  The  correlated  dipole  moment  and  polari¬ 
zability  are  found  to  be  in  excellent  agreement  with 
experiment,  although  the  anisotropy  is  somewhat 
farther  away.  The  hyperpolarizabilities  seem  to 
be  significantly  affected  by  correlation,  and  al¬ 
though  this  is  only  a  single  molecule  out  of  many, 
it  suggests  that  a  theoretical  approach  that  at¬ 
tempts  to  predict  and  explain  the  experimental 
values  for  hyperpolarizabilities  must  definitely 
take  into  account  the  effects  of  correlation. 

In  Table  IV  are  listed  the  individual  correlation 
corrections  for  HF  which  are  of  some  interest  in 
answering  questions  about  the  order  of  perturba¬ 
tion  theory  needed  to  get  converged  correlated 
answers  as  well  as  the  effect  of  the  different  types 
of  diagrams.  It  is  apparent  that  most  of  the  cor¬ 
relation  correction  is  obtained  from  just  the 
second-order  energy  diagrams.  In  fact,  a  com¬ 
parison  of  the  second-order  results  with  the  SDQ- 
MBPT(4)  values  demonstrates  that  there  is  not 
too  much  change  due  to  the  third-  and  fourth-order 
diagrams,  which  are  of  opposite  sign.  In  general, 
the  fourth-order  terms  have  a  somewhat  larger 
magnitude  than  third  order,  further  augmenting 
the  second-order  result.  This  is  particularly  true 
This  behavior  has  also  been  observed  in 
studies  of  molecular  correlation  energies,25  as 
illustrated  currently  by  W(0)  for  the  HF  molecule. 
At  first  sight,  this  may  cause  some  reservations, 
but  comparison25  between  the  DQ-MBPT(4)  model 
for  the  correlation  energy  with  the  infinite-order 
sum  of  double  and  quadruple  excitation  diagrams 
normally  known  as  the  coupled-cluster-doubles 
(CCD)28  model  shows  agreement  to  within  a  milli- 
hartree  for  a  large  group  of  atoms  and  molecules. 
Hence  even  though  fourth  order  is  larger  than 


third  order  for  this  example,  it  d'>es  nit  necessa¬ 
rily  follow  that  higher-order  terms  will  be  im¬ 
portant.  However,  the  very  small  differences 
involved  in  finite-field  calculations  might  be  more 
likely  to  be  affected  by  higher-order  terms:  hence 
this  question  should  be  borne  in  mind.  We  will 
resolve  this  point  in  future  work. 

For  each  property  except  the  energy,  the  single- 
excitation  diagrams  provide  a  larger  fourth-order 
contribution  than  the  double-excitations,  while  the 
quadruple  excitations  are  typically  a  factor  of  2 
to  3  smaller  than  the  fourth-order  double-excita¬ 
tion  contribution.  This  reflects  the  residual  im¬ 
portance  of  single-excitation  terms  for  properties 
dependent  upon  one-electron  ope  rators.  even  though 
this  type  of  single-excitation  contribution  only 
appears  tbecause  of  Bnllouin  s  theoreml  in  the 
fourth-order  energy  via  their  coupling  through 
double  excitations.  The  predominant  single-excita¬ 
tion  effect,  as  discussed  previously,  is  included 
at  the  CHF  level.  Triple-excitation  diagrams 
also  occur  in  fourth  order  and  are  likely  to  be 
somewhat  more  important  than  the  quadruple  ex¬ 
citations.  These  probably  provide  a  correction 
with  the  same  sign  as  the  fourth-order  single-  and 
double-excitation  diagrams.  The  relative  unim¬ 
portance  of  the  quadruple  excitations  plus  the 
fact  that  tae  quadruple  and  triple  excitations  should 
have  an  opposite  effect  on  a  polarizability  suggest 
that  a  limitation  to  single-  and  double-excitation 
diagrams  possibly  summed  to  all  orders  should 
provide  reliable  results  for  these  properties,  if  a 
well-chosen,  well-balanced  basis  set  is  used. 
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ABSTRACT 

Correlated  calculations  using  many-body  perturbation  theory  and 
good  basis  sets  for  the  dipole  moment,  polarizability,  and  hyperpolarizabil ities 
of  t^O  are  reported.  The  finite  field,  coupled  Hartree-Fock  theory  is  used 
as  a  zeroth-order  approximation,  with  electron  correlation  included  via  the 
linked-diagram  theorem.  Single-,  double-,  and  quadruple-excitations  contribu¬ 
tion  are  included.  It  is  found  that  correlation  changes  |b(w=0)|  by  about 
50  percent  compared  to  the  CHF  result.  Even  larger  changes  are  found  among 
the  componts  of  y.  The  dependence  of  the  6  hyperpolarizability  on  geometric 

%  "'ll 

displacements  is  investigated,  finding  that  B  can  be  highly  sensitive  to  slight 
changes  in  bond  lengths. 


I. 


INTRODUCTION 


In  previous  work,^  hereafter  called  I,  a  correlated  study  of 
the  HF  molecule  employing  many-body  perturbation  theory  (MBPT)  of  the  dipole 
moment,  y_,  and  the  polarizability  and  hyperpolarizability  tensors  a,  8,  and  y 
was  reported.  The  objective  of  that  study  was  to  investigate  the  dependence 
of  hyperpolarizabilities  on  normal  basis  sets  of  contracted  gaussian  type 
(CGTO)  and  to  assess  the  effect  of  correlation  on  the  predictions  of  these 
quantities.  Basis  sets  were  chosen  by  using  the  criteria  that  good  agreement 
should  be  achieved  between  the  predictions  of  the  finite  basis  sets  and  numeri¬ 
cal  Hartree-Fock  calculations  of  the  z-components  of  the  higher  polarizabil i- 
f  1  21 

ties.  *  Although  other  components  are  not  sampled,  this  procedure  hopefully 
provides  some  reason  to  believe  that  the  basis  sets  used  would  at  least  be 
adequate  to  provide  a  reasonable  estimate  of  the  correlation  effects  in  hyper¬ 
polarizabilities.  In  I,  correlation  was  found  to  change  pz,  a,  B  =  (&zzz  + 

B  +  B  _)  and  Y  by  6.5%,  14%,  28%  and  39%  respectively,  which  suggests 
that  non-correlated  predictions  of  hyperpolarizabilities  are  unlikely  to  pro¬ 
vide  reliable  agreement  with  experiment. 

In  the  case  of  HF,  no  experimental  values  are  available,  hence  it 
is  of  interest  to  apply  the  techniques  developed  in  I  to  make  predictions  of 
the  hyperpolarizabilities  of  HjO.  Experimental  values  obtained  by  dc  -induced 
second  harmonic  generation  experiments^  are  available  for  H20^  yet 
previous  ^5-7)ion  -correlated  theoretical  calculations  have  not  been  able  to 
provide  reasonable  agreement. 

In  this  paper  we  report  a  series  of  computations  in  different  basis 
sets  for  the  H2O  molecule  ranging  in  size  from  50  to  75  GCTO's.  In  addition, 
we  study  the  dependence  of  hyperpolarizabilities  on  vibrational  stretching  and 


2 


bending  motions,  and  find  a  surprising  degree  of  dependence  of  s  on  the 
symmetric  stretching  mode  of  1^0. 
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II.  SYNOPSIS  OF  THEORETICAL  CALCULATIONS 

The  notation  and  details  of  the  theory  are  reported  in  I.  Briefly, 
we  consider  the  energy  in  the  presence  of  an  electric  field  E_, 

W(E)  =  W(o)  -  -  (2!)  1 “i 3ijk^iEjEk 

"  (^!)  YijklEiEjEkEl  *  *  ‘ 

where  summation  over  repeated  indices  is  assumed.  The  induced  dipole  moment 

is  obtained  as  the  derivative,  -(-=■  ). 

3ti 

In  this  work  finite  field  techniques  are  used,  where  W(£)  is  com¬ 
puted  at  a  series  of  field  strengths  from  which  the  various  components  in 
the  expansion  in  Eqn.  (1)  can  be  obtained.  The  finite  field  strengths  are 
chosen  symmetrically  (e.g.  -  0.01)  which  enables  one  to  eliminate  the  next 
higher  term  in  the  series,  and  thereby  the  principal  contaminating  effect,  by 
taking  odd  and  even  combinations  of  the  finite-field  energies.  H) 

The  energies  are  obtained  from  the  1 inked-diagrams  theorem  of  MBPt(8-10) 

to  be 

oo 

W(E)  -  Wchf(E)  +  l  <  tf0(£)|V(E)[(W0(E)  -  F(E))'1V(E)j<!0o(E)>L  (2) 

K=1 

F(E_),  the  Fock  hamiltonian  plus  the  one-electron  perturbation  due  to  the 
electric  field,  serves  as  the  unperturbed  problem,  while  V(£)  is  the  correlation 
perturbation.  V(£)  is  also  field  dependent  due  to  its  including  the  effective 
SCF  potential.  WQ(£)  is  the  sum  of  the  field  dependent  orbital  energies 
associated  with  the  one-electron  SCF  hamiltonians  that  compose  F(£).  <$0(£) 

is  the  finite-field  SCF  solution.  The  unperturbed  result  in  this  scheme  is 
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the  coupled  Hartree-Fock  (CHF)  solution.  Notice,  in  a  double  perturbation 
approach  relative  to  the  field-free  SCF  solution,  the  CHF  result  corresponds 
to  the  sum  of  an  infinite  series  of  MBPT  diagrams^.  Hence,  the  present 
approach  deals  with  "dressed"  diagrams  relative  to  the  double  perturbation 
approach.  In  this  work,  all  possible  diagrams  that  arise  through  fourth-order 
that  involves  single,  double,  or  quadruple  (i.e.  fourfold)  excitations  are 
included  in  the  calculation  of  W(E^),  while  those  fourth-order  diagrams  that 
have  threefold  excitations  are  neglected.  This  model  is  referred  to  as 
SDQ-MBPT(4) . 0 
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III. _ RESULTS  AND  DISCUSSION 

A.  Previous  Results 

The  electric  susceptibility  for  H90  has  been  determined  from 

11  1 

dc-induced  second-harmonic  generation  experiments  by  Ward  and  Miller  to  be 
( 41 

-94+4  esu/molecule.  '  This  value  is  frequency  dependent,  but  it  is  expected 
that  this  result  should  be  within  about  10  percent  of  the  static  value.  Hence, 

(2) 

x,,  -v  10/3  B(u>-0)  =  10/3(6  +6  +6  ) .  The  minus  sign  signifies  that  8  is 

1  1  zzz  aaz  yyz 

opposite  in  direction  to  the  permanent  dipole  moment  of  H^O,  assumed  to  be  in 
the  direction  O'H^  . 

Several  previous  calculations  of  6  at  the  uncoupled  Hartree-Fock  (UCHF) 
and  coupled  Hartree-Fock  (CHF)  level  have  been  reported. ^5'7^  Liebman  and  Moskovitz 
obtain  +90.6  in  a  UCHF  calculation  for  where  the  sign  seems  to 

be  in  error.  Arrighini,  et.  al .  ^  report  CHF  values  of  -51.6  and  -48.0 

in  (2sl p 1 1 s )  and  (5s3pld|2slp)  Slater  (STO)  basis  sets,  respectively.  In  the 
larger  STO  basis,  the  hydrogen  components  perpendicular  to  the  plane  are  omitted. 

It  is  a  common  result  that  |sj  is  smaller  when  a  better  basis  set  is  used,  although 
the  reduction  in  |b|  is  usually  more  dramatic^’7)  than  is  this  example. 

Lazzaretti  and  Zanasi  report  CHF  calculations  of  6  for  H£0  for  three 
contracted  gaussian  orbital  (CGTO)  basis  sets!7  ^  ranging  from  minimum  (2slp|ls), 
to  double  zeta  (4s  2p|2s),  to  a  polarized  double  zeta  (4sl  2pl  d|2sl  p)basis  set. 

The  values  for  xj^(w=o)are  respectively,  -52.5,  -79.2,  and  -21.9.  Again  the 
general  trend  that  a  better  basis  reduces  CHF  values  of  |s|  is  evident,  emphasiz¬ 
ing  the  accidental  agreement  between  experiment  and  the  double-zeta  basis. 

(it) 

Unfortunately,  the  very  accurate,  large  basis  results  of  Werner  and  Meyer 


6 


for  the  dipole  moment  and  polarizability  are  not  extended  to  predictions  of 
hyperpolarizabilities,  since  these  authors  correctly  question  whether  even 
basis  sets  much  larger  than  those  used  in  the  previous  studies  can  be  expected 

(  1  o  ) 

to  have  any  validity  for  properties  as  sensitive  as  hyperpolarizabilities. 

In  I  we  showed  that  correlation  and  basis  sets  can  have  a  very 
large  effect  on  ab  initio  predictions  of  hyperpolarizabilities,  so  it  is 
of  interest  to  use  the  same  techniques  to  attempt  to  make  some  assessment 
of  the  capability  of  theory  to  accurately  predict  hyperpolarizabilities. 

The  present  study  will  focus  on  three  elements:  correlation,  basis  sets, 
and  the  dependence  of  6  on  changes  in  molecular  geometry. 

B.  Basis  Sets 

In  this  work,  four  different  basis  sets  are  considered,  all  larger 
than  those  previously  used  in  studies  of  H20  hyperpolarizabilities.  The 

) 

prescription  used  is  essentially  that  of  Christiansen  and  McCullough. 

The  first  three  basis  sets,  A,  B,  and  C  are  built  from  Dunning's  5s3p  contraction^ 

( i  c ) 

of  Huzinaga's  primitive  9s5p  basis  set  for  the  oxygen  atom;  'and  Dunning's  3s 
contraction  for  hydrogen.  These  functions  are  augmented  by  a  3d  gaussian 
polarization  function  on  oxygen  (all  six  components)  (exponent  -  1.211)  and 
a  p  gaussian  polarization  function  (exponent  -  0.761)  on  hydrogen.  This 
(5s3pld(3slp)  basis  is  considered  to  be  generally  adequate  to  describe  the  basic, 
molecular  charge  density,  but  insufficient  for  predictions  of  polarizabilities. 

To  assist  in  describing  such  additional  properties,  the  basis  is  augmented  by 
adding  an  extra  set  of  s,  p,  and  d  functions  to  oxygen  with  the 
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exponents  chosen  to  have  the  same  ratio  as  the  last  two  functions  of  each 
type,  generally  a  factor  of  about  0.3  times  the  previous  exponent.  Similarly, 
hydrogen  is  augmented  by  an  s  and  a  p  function  whose  exponents  are  chosen  in 
the  same  manner.  This  results  in  basis  A,  a  (6s4p2d |4s2p)  basis  of  50  CGTO. 

It  is  expected  that  this  basis  should  be  generally  reliable  for  polarizabilities, 
buy  not  hyperpolarizabil ities  where  at  least  another  shell  of  functions  with 
increasingly  diffuse  exponents  is  recommended.  By  proceeding  in  this  manner, 
an  additional  set  of  s,p,  and  d  functions  are  added  to  oxygen,  and  an  additional 
s  and  p  function  to  hydrogen.  This  constitutes  basis  B,  a  (7s5p3d|5s3p)  basis 
of  68  CGTO.  To  add  additional  diffuseness  to  the  charge  cloud,  a  third  basis, 
basis  C,  includes  an  additional  s  (exponent  0.0079)  and  p  function  (exponent 
0.0057)  located  at  the  center-of-mass  of  the  molecule.  This  (7s5p3d|5s3p;lslp) 
basis  requires  72  CGTO. 

The  construction  of  the  basis  sets  A,  B,  and  C  emphasizes  the  long- 
range,  tail  region  of  the  molecule  where  more  and  more  diffuse  functions  are 
employed  to  try  to  describe  this  part  of  the  charge  density.  In  I  we  found 
that  by  comparison  with  numerical  Hartree-Fock  calculations  of  the  on-axis 
components  of  the  polarizability  and  hyperpolarizabilities,  it  is  also  important 
to  have  a  basis  that  is  capable  of  describing  the  small  differences  in  the 
charge  density  in  regions  moderately  near  the  nuclei.  Since  our  basis  in  I 
is  successful  in  obtaining  quite  good  agreement  with  the  numerical  Hartree-Fock 
results,  we  also  have  used  this  procedure  to  construct  a  fourth-basis  set,  D, 
that  we  hope  will  retain  some  of  the  apparent  reliability  of  the  previous  work. 

This  basis  starts  with  Dunning's  (4s3pj2s)  contraction,  then  uses  a 

(14) 

d-polarization  function  defined  as  a  two  gaussian  fit'  '  to  a  Slater  exponent  of 
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3.0175.  The  smaller  exponent  is  then  multiplied  by  0.4  to  generate  a  new 

d-orbital  exponent.  Two  successive  sets  of  d-functions  are  obtained  by  a 

repetition  of  this  scheme.  Similarly,  a  fourth  p-function  (exponent  0.2137) 

is  added  to  oxygen,  and  two  more  p-functions  with  exponents  chosen  to  be 

0.4  times  and  then  (0 . 4 )2  times  this  value.  Also,  two  s-functions  are  added 

by  taking  0.4  and  (0.4)2  times  the  smallest  s-exponent  (i.e.,  0.2846).  Two 

additional  s-functions  are  added  to  hydrogen  (exponents  0.1776  and  0.04932) 

(14) 

and  a  2p  polarization  function  two  gaussian  fit'  to  an  ST0  exponent  of  2.047> 
and  an  additional  p-function  with  exponent  (0.2174).  The  resulting  basis  of 
65  CGT0  is  of  the  form  (6s5p4d  1 4s2p ) .  The  differences  in  the  construction  of 
this  basis  compared  to  A,  B,  and  C  should  not  be  major,  but  are  of  interest  in 
assessing  the  sensitivity  of  the  predictions  to  modest  changes  in  basis  set, 
and  for  comparisons  with  the  previous  results  for  HF.^ 
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C.  Results 


In  Table  I  are  listed  CHF  results  for  u_,  3,  and  %  for  the  basis 

sets  A-D.  Since  bases  B  and  C  are  constructed  from  basis  A  by  simply 

adding  additional  diffuse  functions  without  any  other  change,  it  is  possible 
to  form  some  assessment  of  the  effect  this  augmentation  has  on  the  predicted 
results. 

The  largest  change  is  found  in  the  components  of  !b| , reducing  their 
size  as  has  been  observed  to  usually  be  the  case  when  a  better  basis  set  is 
employed.  Also,  since  the  secona-harmonic  generation  experiment  gives  a 
| 0 |  larger  than  the  theoretical  predictions,  improving  the  basis  set  further 
destroys  any  fortuitous  agreement  between  the  CHF  results  and  experiment. 

The  highly  diffuse  functions  included  in  basis  C,  but  not  B,  are  observed  to 
have  only  a  small  effect  on  £  . 

A  significant,  but  smaller  change  occurs  in  the  values  of  a  between 
basis  A  and  B,  although  the  dipole  moment  is  little  modified.  Again,  the 
extra  diffuse  functions  in  C  have  little  effect  on  a  and  none  on  y.  The 
change  in  a  between  A  and  B  is  consistent  with  the  observation  of  Werner  and 
Meyer^3^(WM)  that  at  least  a  3d  basis  set  is  necessary  to  simultaneously  describe 
the  energy,  dipole  moment,  and  polarizability.  In  fact,  for  HgO  WM  recommend 
the  d-orbital  exponents  0.10,  0.30,  and  1.20,  while  our  basis  B  uses  the 
quite  similar  though  independently  arrived  at  values  of  0.11,  0.37,  and  1.211. 

The  main  problem  in  basis  sets  A,  B,  and  C  is  found  in  the  dipole 
moment.  Although  there  is  little  change  when  more  diffuse  functions  are 
added,  its  value  is  still  rather  far  from  the  WM  CHF  result  of  uz  =  0.782. 
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Consequently,  the  failure  of  bases  A-C  to  account  for  this  resi.lt,  appears 
to  lie  among  the  interior  basis  functions  rather  than  the  more  diffuse 

functions.  Since  B  tends  to  be  sensitive  to  both  interior  and  tail  regions 

Or 

of  the  charge  density,  in  hopes  of  obtaining  more  reliable  results  for  _s 

we  focus  on  a  slightly  modified  approach  to  generating  a  basis  set  for 

hyperpolarizabilities  that  we  found  to  be  successful  for  HF,^  which 
gives  basis  set  0.  This  basis  uses  a  two  gaussian  fit  for  the  tight 
d-function,  and  an  additional  diffuse  d-orbital,  as  well  as  some  other, 
presumably  less  significant  changes.  In  Table  I  it  may  be  seen  that  for 

basis  D,  y,  and  a  are  in  excellent  agreement  with  the  WM  results  of 

'V  % 

y  =  0.782,  and  a  =  8.47,  a  =  9.04,  and  a  =  7.99.  As  observed  by  WM 

XX 

the  polarizability  is  less  sensitive  to  such  a  change  than  is  the  dipole  moment. 

Again,  as  in  the  transition  from  A  to  B  and  C,  basis  D  further 
reduces  the  components  of  6  .  It  is  not  appropriate  to  claim  any  basis  set 

*\j 

convergence  from  these  four  sets  of  CHF  calculations ,  but  experience  strongly 
indicates  that  some  basis  set  limit  convergence  would  tend  to  predict  even 
smaller  jej  at  the  CHF  level,  further  destroying  any  agreement  between  the 
CHF  predictions  and  experiment.  To  rectify  this  dilemma,  three  features  can 
be  addressed  that  can  increase  |b|.  These  are  correlation,  geometric  displace¬ 
ment  as  might  become  important  in  higher  temperature  experiments,  and  the 
frequency  dependence.  In  this  paper  we  will  consider  the  first  two  possibilities 
with  frequency  dependence  to  be  studied  in  future  work. 

Considering  the  second  possibility,  first,  Figures  1  and  2  illustrate 

the  dependence  of  two  components  of  B  on  the  symmetric  stretch  and  bending 

2  1/2 

mode  in  HgO.  These  CHF  calculations  use  basis  C.  The  values  <aRj>  '  » 


0.128b 
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1/2 

<a02>  =  8.72  are  the  experimental  root-mean-square  amplitudes  for  H^O 

displacement  coordinates,  at  °K. The  values  of  <aR-j >  =  0.0266  and 

<A0>  =  0.183  degrees.  It  is  apparent  that  there  is  a  very  steep  dependence 

of  6  and  6  on  the  stretch,  although  not  the  bending  mode.  In  the  case 

of  Bxxz»  its  value  is  comparatively  unaffected,  changing  from  +0.26  at 

r  +0.128  a.u.  to  -0.68  at  r  and  +1.15  at  r  -0.128  a.u. 
e  e  e 

Judging  from  Figures  1  and  2,  the  zero-point  correction  of  b  1 
(3  +B  +3  )  amounts  to  about  1.7  units.  The  sensitivity  of  B  and 

£  AAA  V  V  4u  L.  Z. 

6^^  to  displacements,  however,  suggests  that  experiments  conducted  at 
elevelated  temperatures  involving  higher  rotational  states,  should  be 
correlated  for  excited  state  populations  and  centrifugal  distortion  of 
bond  lengths. 

The  change  in  e  with  displacement  is  more  extremem  than  is  the 

change  in  a.  The  mean  polarizability  assumes  the  values  8.54  at  r  ,  and  7.76 

at  R  -oR,  and  9.47  at  R  +6R  .  In  the  case  of  the  bending  mode,  a  =  8.52  at 
e  e  - 

0  -60  and  8.60  at  0  ^9. 
e  e 

The  other  element  that  is  important  in  making  accurate  predictions 
of  static  hyperpolarizabilities  is  electron  correlation.  In  Table  II  are 
listed  the  CHF  and  correlated  SDQ-MBPT(4)  results  for  basis  sets.  A  and  D, 
with  experimental  values  where  available. 

The  effect  of  correlation  changes  the  CHF  dipole  moment  in  basis  D 
by  -^6  percent,  a  by  ^9  percent  and  g  by  ^50  percent.  Similarly  large  changes 
in  the  various  components  of  %  of  -v36  to  85  percent  are  observed.  In  the 
case  of  HF,  similar  values  of  ^6  percent,  ^14  percent,  28  percent,  and  39-to- 
75  percent  are  found.  ^  In  v  ,oii  6,  and  y  the  effect  of  correlation  is 
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approximately 
case  of  u,  a, 

—  "V 

additive,  but 


the  same  as  the  basis  set  effect  between  bases  A  and  D.  In  the 
and  y  the  basis  set  effect  and  the  effect  of  correlation  are 
for  8  this  is  not  at  all  the  case. 


The  delicate  balance  between  basis  set  and  correlation  needed  to 
describe  6  accurately,  has  mixed  effects.  Since  a  larger  correlation 
correction  will  typically  be  obtained  in  a  better  basis  set,  as  the  basis  is 
improved,  the  observed  reduction  in  the  CHF  | s |  will  probably  be  largely 
offset  by  the  greater  correlation  contribution.  This  has  the  implication 
that  some  hope  exists  for  obtaining  a  reasonable  level  of  convergence  with 
basis  sets  for  predictions  of  |e|.  However,  at  the  CHF  level,  no  such 
possibility  persists.  On  the  other  hand,  the  sensitivity  of  8  to  both  basis 
sets  and  correlation  emphasizes  the  importance  of  treating  both  aspects  of 
the  problem  in  a  balance  fashion,  and  at  a  higher  level  of  sophistication. 

The  results  of  basis  D  in  Table  II  correspond  to  an  electric 

susceptibility  xj^)U=0)  of  -66  at  rg  and  -74  at  rQ,  still  about  21  percent 

different  than  the  frequency  dependent  SHG  results  of  -94+4.  It  is  clear, 

however,  that  what  agreement  there  is  between  experiment  and  prior  CHF 
(5-7) 

calculations  '  has  to  be  essentially  accidental. 

In  Table  III  are  reported  the  individual  correlation  corrections  to 
£,  8,  and  y.  As  in  the  prior  study  of  HF,  second-order  perturbation  theory 
provides  most  of  the  correlation  contribution.  Since  in  the  finite-field 
approach  used  in  this  work,  most  of  the  single  excitations  effects  that  are 
quite  important  for  one-electron  operators  are  incorporated  into  the  CHF 
(  unperturbed)  results,  most  important  correlation  effects  are  due  to 
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double-excitation  relative  to  field-dependent  orbitals,  which  are  exclusively 
responsible  for  second-  and  third-order  perturbation  theory.  In  fourth- 
order,  the  residual  effects  of  single  excitations  first  appear,  and  their 
value  is  slightly  larger  than  the  fourth-order  double  excitation  contribution. 
As  in  I,  the  quadruple  excitation  fourth-order  diagrams  contribute  very  little 
to  the  correlated  result. 

In  summary,  the  importance  of  electron  correlation  to  predictions 
of  hyperpolarizabilities  has  been  emphasized.  Correlation  effects  make 
a  contribution  that  significantly  increases  |e|  and  the  components  of  y. 

'V 

Since  the  experimental  values  for  jf$|  are  normally  larger  than  the  theoretical 
predictions,  correlation  can  make  an  important  contribution  toward  resolving 
this  discrepancy.  It  is  fairly  clear,  that  CHF  results  can  at  best  offer 
fortuitous  agreement  with  experiment  when  inferior  basis  sets  are  used.  We 
have  also  shown  that  the  6  hyperpolarizability  can  be  very  sensitive  to  small 
displacements  of  the  nuclei  in  a  molecule.  This  sensitivity  should  be 
considered  in  future  theoretical  calculations  and  can  be  a  source  of  error 
in  comparing  with  the  experimental  results.  The  remaining  major  question 
to  resolve,  is  the  change  in  s  with  frequency.  That  feature,  too,  should 
increase  the  |e|  obtained  in  theoretical  calculations.  When  this  element 
is  included,  perhaps  reasonable  agreement  between  ab  initio  calculations 
and  experiment  will  finally  become  possible. 
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FIGURE  2.  DEPENDENCE  OF  THE  B22Z  AND  6yy2  COMPONENTS  OF  THE  HYPERPOLARIZABILITIES 
ON  THE  BENDING  MODE  OF  H„0.  CHF  CALCULATIONS,  BASIS  C.  , 


17 


TABLE  I.  DIPOLE  MOMENT  AND  FIRST,  SECOND,  AND  THIRD  ORDER  POLARIZABILITIES 
FOR  H20  DETERMINED  BY  COUPLED-HARTREE-FOCK  THEORY  (Molecule  is  in 
zy-plane,  with  z  the  molecular  axis.  Orientation  is  0*h£  defined  as 
a  positive  dipole  moment;  rQH  =  1.811b,  0  =  104.5°) 


Basis  Set 

A3 

Bb 

Cc 

Dd 

pz 

0.  807 

0.796 

0.796 

0.784 

azz 

8.01 

8.50 

8.55 

8.47 

ayy 

8.78 

9.11 

9.12 

9.18 

axx 

7.21 

7.96 

7.95 

7.95 

ezzz 

-8.08 

-6.65 

-6.31 

-5.22 

syyz 

-12.86 

-10.03 

-9.75 

-9.63 

r 

0. 

0.85 

-0.68 

-0.48 

^xxz 

Y 

490 

_ 

770 

zzzz 

V 

_ 

_ 

270 

■yyzz 

V 

• 

_ 

350 

xxzz 

Y 

_ 

1400 

xxxx 

Y 

_ 

_ 

460 

yyyy 

_ 

'xxyy 

a  Basis  A  is 
b  Basis  B  is 
c  Basis  C  is 
d  Basis  D  is 


(6s4p2d |4s2p) ,  =-76.05391. 

(7s5p3d |5s3 p) ,  =  076.05412. 

(7s5p3d|5s3p;lslp);  =  -76.05420. 

(6s5p4d |4s2p) ;  =  -76.05443. 
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Footnotes  (Table  II.) 

a  Calculations  in  basis  sets  A  and  0  are  at  r0  =  1.811  a.u.  and  0  =  104.5 


